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■Abstract 

In this work we investigate the properties of the broad emission line components in the Balmer series of a sample of 
.Type 1 Active Galactic Nuclei ( AGN) . Using the Boltzmann Plot method as a diagnostic tool for physical conditions in 
the plasma, we detect a relationship among the kinematical and thermo-dynamical properties of these objects. In order 
to further clarify the influence of the central engines on the surrounding material, we look for signatures of structure in 
the broad line emitting regions, that could affect the optical domain of the observed spectra. Using a combination of line 
.profile analysis and kinematical modeling of the emitting plasma, we study how the emission line broadening functions 
are influenced by different structural configurations. The observed profiles are consistent with flattened structures seen at 
quite low inclinations, typically with i < 20°. Since this result is in good agreement with some independent observations 
at radio frequencies, we apply a new formalism to study the properties of AGN central engines. 
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1. Introduction 

The spectrum of a Type 1 AGN is characterized by 
.the presence of prominent broad emission lines, with a full 
"width at half the maximum (FWHM) which corresponds 
■to velocity fields exceeding lO'^kms^^. The properties 
of the emission lines are connected to the physical con- 
ditions within the source and they suggest that the lines 
are mainly originated by electr on - ion recombinat ions in a 
photo- ionized plasma (see e. g. lOsterbrock . TqSQI ). Unfor- 
tunately, since the broad components of these spectral fea- 
tures are originated very close to the center of the source, 
the structure of the so called Broad Line Region (BLR) 
cannot be spatially resolved by our observations. This 
'limit affects our ability to investigate the physics of AGN 
and to provide a self-consistent theoretical interpretation 
of their observational properties. 

At present it is widely accepted that AGN are powered 
by matter accreting into the gravitational field of a Super 
Massive Black Hole (SMBH), since the gravitational bind- 
ing energy of the accretion flow may provide the ionizing 
radiation fleld, required to drive the observed line emis- 
sion from the surrounding plasma. Assuming the BLR to 
be in virial equilibrium under the d ynamical influence of 
this central engine (see, for instance. iPeterson fc Wandel 
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19991: IWoo fc Urrv I . [20o3 ISulentic et al. I . l2006l etc.) the 
mass of the black hole might be estimated with an expres- 
sion in the form of: 



Mbh = / 



G 



(1) 



where G is the gravitational constant, -Rblr an estimate 
of the BLR extension, and Av a measurement of the gas 
radial velocity distribution. The factor / accounts for the 
unknown structure of the gas motion pattern. It is re- 
quired in order to estimate the actual kinematical prop- 
erties from the line of sight velocity distribution, that we 
measure in the spectra, and it is therefore referred to as 
the geometrical factor. 

Many works deal with the relationship among the emis- 
sion line pro files and the actual dynamics of the BLR (see . 



198nl . Il98ll : iNetzer 1 . 119901 : 



for insta nce. iCaoriotti et al. 

iFerland et al. I . ll992h Recen t in vestigations, such as thos e 



of IVestergaard et al. I (|2000[) and lNikolaiuk et al. I (|2005l ). 

point towards a considerably flattened geometry, which is 
consistent with the existence of an accretion disk around 
the SMBH. In this progress report, we describe the analy- 
sis which we performed on a sample of Type 1 AGN spec- 
tra, in order to examine the influence of the AGN central 
engines onto the surrounding BLR plasma. With a combi- 
nation of observations, analytical calculations and theoret- 
ical models, we look for indications of the BLR structural 
properties in the emission line profiles. We organize this 
report as follows: in §2 we introduce the relation among 
BLR kinematical and thermo-dynamical properties; in §3 
we outline the analytical formalism, that we exploit to 
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Figure 1: Four examples of Boltzmann Plots applied to the spectra of 
the BLR in some Type 1 AGN: in the upper left panel no straight line 
fit is achieved (class iv); in the upper right panel only a poor fit, with 
a relative uncertainty AA/A > 0.2 in the temperature parameter, is 
performed (class iii); the bottom left panel gives a good fit, but He 
could not be detected in the spectrum (class ii); finally the bottom 
right panel shows a straight line fit to the normalized intensities of 
the Balmer series up to He (class i). 



investigate the BLR structure; in §4 we discuss the predic- 
tions of our BLR structural models; finally §5 is devoted to 
the discussion of the results achieved by comparing mod- 
els and observations, with a mention to some possible tests 
and to the implications related to the use of spectroscopic 
measurements in our sample of AGN. 

2. Plasma physics with the Boltzmann Plots 

Due to its high density and strong interactions with the 
ionizing radiation field, the BLR plasma cannot be stud- 
ied by means of the standard spectroscopic techniques, 
that are commonly adopted in other low density nebular 
environments. This characteristic forces us to adopt al- 
ternative techniques to investigate its properties. Looking 
at the flux in the broad component of an emission line, it 
is possible to introduce a normalized line intensity^ with 
respect to the atomic constants which characterize the cor- 
responding transition, in the form of: 



Aui Fui 
9uA 



ul 



(2) 



where A„; is the line wavelength for a transition from an 
upper level u to a lower level I, Fui is the observed flux, A^i 
the spontaneous transition probability, and gu the upper 
level's statistical weight. Looking at the normalized emis- 
sion line intensities of some spectral features, belonging 



Figure 2: A parameter as a function of the BLR velocity fields, 
inferred by both the FWHM and FWZI of the Balmer lines: the 
open circles are sources of class i; the open squares represent objects 
of class ii; filled triangles are spectra of class iii and, finally, asterisks 
are objects with BP class iv. The crosses in the lower right corners 
of each panel are the median uncertainties of measurements. 



to the same transition series, it can be shown that if the 
plasma approaches a conditi on of Local Thermo-dy namical 



Equihbrium (LTE), we get (jPopovic I . l2003l . l200i 



logIn = B-AEu, 



(3) 



where Eu represents the upper level's excitation energy for 
the transition, B is a constant for all the emission lines in 
the series, while A takes the role of a temperature param- 
eter, corresponding to: 



A = loge/(fcBTe), 



(4) 



with fee representing the Boltzmann constant and T^. an 
average value of the electron temperature in the plasma. 

The diagrams which show the emission line normalized 
intensities of the Balmer series as a function of the upper 
level's excitation energies are known as Boltzmann Plots 
(BP). A strict interpretation of the BP results, however, 
requires that the physical conditions in the line emitting 
region should not change abruptly. Otherwise the sim- 
ple formulation of the emission line intensities in the se- 
ries would not be possible and the LTE approximation 
would subsequently fail. In Fig. 1 we give some examples, 
showing four different situations, that we commonly ob- 
serve, applying the BP to the BLR plasma of some AGN 



( La Mura et al. 1 . 120071 ). Extracting the BLR signal from 



our spectra, we computed the emission line normalized in- 
tensities of the Balmer series as in Eq. (2). We performed 



2 



several flux measurements, that were subsequently aver- 
aged together and compared with the spectral noise fluc- 
tuations, in order to estimate the mean values and their 
uncertainties, therefore assuming the error of the BP to be 
given by: 

A(log/„) = ^, (5) 

with apui representing the emission line flux dispersion. 
Repeating the emission line flux estimates was a funda- 
mental step, because of the large difficulties associated 
with flux measurements in the broad line components, in a 
spectral region where several contributions to the observed 
emission are blended together. 

The temperature parameters extracted from our BP 
and illustrated in Fig. 2 appear to be related to the emis- 
sion line profiles, in the sense that higher values of A and 
better realizations of the prediction of Eq. (3), are com- 
monly found in the range of broad line emitting objects. 
Unfortunately our sample is not particularly representa- 
tive in this region, but the circumstance suggests that the 
thermo-dynamical properties of sources with various emis- 
sion line widths arc under the infiuence of different ionizing 
radiation fields, as it would be expected when powering the 
central engines of these AGN with different combinations 
of SMBH mass and accretion rate. In order to improve our 
understanding of the role played by these fundamental pa- 
rameters, we further investigated the processes which lead 
to the formation of the observed line profiles and their 
relations with the actual source's dynamical properties. 

3. Balmer emission line profile analysis 

3.1. Emission line broadening from cross- correlation 

The BLR spectrum shows several emission lines cor- 
responding to many permitted and some semi-forbidden 
transitions of variously ionized atomic species. In the op- 
tical domain, the Balmer series of hydrogen provides a set 
of bright emission lines belonging to a well defined family 
of interactions among matter and radiation, whose profiles 
are infiuenced by the kinematics of gas. 

To perform our study, we adopt ed the cross-correlation 
formalism, originally described by Tonrv fc Davis ( 1979f l 
and then updated bv lStatler" ( 1995| ). As a starting point, 
we can approximate the spectrum of a broad line emitting 
source as the convolution of an appropriate template of 
narrow emission lines Tix) with the corresponding broad- 
ening function B{x): 



S(x) ~T[x)*B(x), 



(6) 



where S{x) is the observed spectrum, while x represents a 
logarithmic wavelength coordinate of the form x — Ahi A-l- 
B, such that the effect of radial velocities results in linear 
shifts along x. Introducing the cross-correlation function 
of the spectrum with the template: 

X{x) ^ S{x)®T{x) ^ [ S{x)T{x + x')dx\ (7) 



it can be shown that the cross-correlation function (XCF) 
approximates the convolution among the template's auto- 
correlation function (ACF) an d the kinemat i cal broaden- 



ing fu nction (BF) of the object (jStatler 

X{x) ~ [T{x)®T{x)] * B{x). 



1995: La Mura et al. 



(8) 



Since T{x) is known and X{x) is drawn from observations, 
as far as the template is correct, it is possible to recover 

Restricting our analysis to the primary cross-correlation 
peak, which carries most of the kinematical information 
and it is weakly affected by template mismatch, Eq. (8) 
has a discrete form that, using the simplified notation 
Fi = F{xi), becomes: 



Xk 



N 

E 



N 



Bk 



(9) 



Provided that all the functions are null when they are com- 
puted outside the range < i < A^, Eq. (9) defines a 
system of A^ + 1 linear equations in the A^ -|- 1 variables 
Bk~i {k> i). A standard x^-minimization routine can be 
therefore used to solve the problem. 

3.2. Analytical expressions for the broadening functions 

The BLR emission lines are influenced by the effects 
of complex kinematics within the source and of radiation 
transfer from the source to the observer. For this reason 
it is hardly conceivable that a simple analytic expression 
might be used to fit the resulting profiles. In the simple 
case of a random motion distribution, the emission line 
profiles would approximately match a Gaussian function, 
but, in presence of ordered kinematical components, we 
expect significant deviations from this shape. A good way 
to estimate the importance of such various effects is to pa- 
rameterize the observed BF by means of a Gauss- Hermite 
orthonormal expansion, similar ly to what is described in 
Van Per Marel fc Franx I (ll993l) for the case of stellar kine- 
matics in elliptical galaxies. Following their method, if we 
call a{v) the normal Gaussian function: 



a{v) = 



1 



cxp 



V 

2^ 



(10) 



where is the line of sight velocity dispersion, we may 
express the emission line BF as: 



B{v) = Boa{v - Vsys) 



N 

4=3 



, (11) 



where we call Bq the BF normalization factor, Vsys the 
systemic radial velocity offset between the BF and the 
chosen reference frame, Hi{v — Vsys) the i*^ order Her- 
mite polynomial, and hi the corresponding coefficient. A 
wide description of the properties of the Hermite polyno- 
mials is given in IVan Per Marel fc Fran^l (|l993l ). It is 
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Figure 3: Mass - inclination degeneracy. In this plot we illustrate 
the expected FWHM in the BF of emission lines originated in a disk 
structure. The empty triangles connected by the continuous line are 
the predicted values for reference disk models with i = 10° in the 
gravitational field of black holes with Mgjj = 1, 2, 4, 6, 8, 10, 12, 
and 14 ■IO'^Mq (respectively models from 1 to 8); the circles with the 
dashed line show the predicted behavior for the reference disk model 
in the gravitational field of the black hole with Mgjj = IO^Mq, when 
seen under inclinations of i = 10°, 15°, 20°, 25°, 30°, 35°, 40°, and 
45° (again from model 1 to 8). 

demonstrated that odd order functions account for asym- 
metric distortions of the Gaussian profile, while even order 
functions have a symmetric effect. Truncating Eq. (11) to 
A'^ = 4, the Hermite polynomials are expressed by: 

iJ3(2/) = -i={2V2y'' - 3V2y) (12a) 

7?4(2/)-^(V-12y2 + 3). (126) 
V24 

Therefore, it is possible to estimate the role of non-Gaussian 
kinematical components, using the whole BF profile, sim- 
ply by fitting the observed shape with a truncated Gauss- 
Hermite series and measuring the appropriate values of hs 
and 



Figure 4: The effect of inclination on the symmetric distortion of 
the BF. Here we show a comparison of the reference model predic- 
tions concerning the BF symmetric component (continuous lines) 
with four variants: in the upper panel we plot models with stronger 
(long dashed line) and weaker (short dashed line) disk emission with 
respect to the surrounding gas distribution; in the bottom panel we 
show the differences obtained by setting o-g^u = 0.07 c (short dashed 
line) and ug^n = 0.09 c (long dashed line). 

a Gaussian function with FWHM = 167 km s""'^. At the 
spectral resolution of Sloan data, the logarithmic sampling 
of the wavelength coordinate can be performed with dis- 
crete bins corresponding to 69 kms~^ each. We compute 
the template's ACF: 

A{x) = T{x) ® T{x) (13) 

and the cross-correlation functions of the BLR spectra 
with the template, following the definition of Eq. (7). Ap- 
plying the least squares formalism to the equation sys- 
tem (9), it follows that the BF of each spectrum must 
satisfy the relations: 

N / N \ ^ 

Y,B^ 51 ^J^^-'^ - ^ (14) 

1=0 \j=0 j i=k 



3.3. The Balmer line broadening functions 

Applying our formalism to a set of BLR spectra, col- 4. BLR kinematical and structural models 
lected at the Sloan Digital Sky Survey (SPSS) spectro- 

scopic database fsee lAdelman-McCarthv et al. l . l2008l : lLa Mura et^ e the BLR structure cannot be represented by a 
20091 for more discussion), our task is then to recover their random motion pattern, the shape of the broad emission 

lines may exhibit large deviations from the Gaussian pro- 
file. If the BLR has a flattened component which is seen 
at low inclination the geometrical structure can consider- 
ably affect the dynamical interpretation of spectroscopic 



BF. To calculate the cross-correlation functions, we build 
a template of Balmer emissio n lines, following themedian 
line intensity ratios found bv lLa Mura et al~T ( 2007 ). The 

template assumes that the SDSS instrumental profile is 
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Figure 5: BF symmetric distortion plotted as a function of 
FWHMh/j. We observe a large evolution of the symmetric com- 
ponent as a function of broadening, suggesting the possibility of in- 
clination effects in nearly face-on flattened structures. Such effects 
become weaker as the profile width increases. 



data. We illustrate this concept in Fig. 3, where we plot 
the expected FWHM in the broadening function produced 
by disks surrounding black holes of increasing mass and we 
compare it to the situation of a black hole of fixed mass, 
but with the disk seen under different inclinations. 

Exploiting a combined BLR structural model, involv- 
ing an accretion disk and a surrounding distribution of gas 



ing an accretion aisK: ana a surrounamg gistribution ot gas 
( Chen et al .111 9891 : Ichen k Halperii1 . ll989l : IPopovic et al. 
2004( 1 ■ we computed a range of emission line profiles, study- 
ing the effect of source inclination. Our reference model as- 
sumes Rin = 1834 i?5 for inner radius, Rblr = 18340 i?5 
for outer radius, aoisk = 0-003 c as the intrinsic velocity 
dispersion in the disk, and a = —2.0 for the radial emission 
power law, where Rs and c represent the Schwarzschild ra- 
dius and the speed of light. Moreover, this model includes 
a surrounding gas distribution having a velocity dispersion 
of aBeii = 0.008 c. This model carries out the best match 
to the observed line profiles with the assumption of various 
disk inclinations. 

In Fig. 4, we compare the reference model with some 
variants, obtained with slightly different parameters. We 
note that all the models predict a strong dependence of the 
symmetric non-Gaussian component (the coefficient /14 in 
the Gauss-Hermite expansion) on the disk inclination, in 
the range of small values of i. The reason is quite simple, 
because a nearly face-on disk enhances the low radial ve- 
locity peak of the BF, while an edge-on disk is more likely 
to affect the high velocity wings. 

As we show in Fig. 5, where we plot the measured 



Figure 6: BLR inclination inferred by comparison among the BF 
symmetric component /14 , predicted by the reference model for vari- 
ous inclinations, and the corresponding distribution observed in our 
data. According to the model predictions, the h4 coefficient is very 
sensitive to fiattening and inclination in nearly face-on structures, 
while it becomes a weaker indicator for larger inclinations. 



values of /14 as a function of FWHMh/j, there is a remark- 
able evolution of the symmetric components of the line 
profiles, which decrease in importance for increasing line 
profile width. Such an effect suggests that a considerable 
variation of the geometrical factor / might be present and 
it should be taken into account in order to estimate the ac- 
tual properties of the SMBH. Using the model predictions, 
we can exploit the symmetric distortion coefficient to esti- 
mate the inclination of the flattened BLR component and 
to apply a correction to our dynamical interpretation of 
the observed line profiles. 

5. Results and discussion 

5.1. The BLR geometry 

It can be shown that completely neglecting the role 
played by the BLR geometrical factor may lead to incor- 
rect black hole mass estimates, with uncertainties that, in 
the worst cases, could span over two orders of magnitude. 
Using the line profile analysis to estimate the BLR struc- 
tural properties, we apply a modified scheme to evaluate 
the black hole mass in Eq. (1), introducing an equivalent 
velocity field, such that Veq — / Aw, defined as: 
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Figure 7: H/3 asymmetry parameter distribution vs. FWHMu/j- 
Filled symbols represent objects where the asymmetric component 
exceeds 10% of the Gaussian contribution. We plot as triangles ob- 
jects that are affected by negative asymmetry, yielding red shifted 
peaks and blue shifted wings, while we use squares to represent pos- 
itive asymmetry sources, with a blue shifted peak and red shifted 
wings. Large asymmetries characterize objects where fits with the 
reference model are more likely to be problematic. The bar in the 
upper left region of the diagram is a median estimate of the mea- 
surement errors. 



Assuming that the line profile broadenin g results from 
both planar and non-planar motions (.Labita et al. 
McLure et a"l~l . l2002t Ijarvis fc McLure I . l2006l etc.). Veq 
combines the velocity estimates obtained from the H/? emis- 
sion line profile by fitting two Gaussian functions, which 
are subsequently compared with the reference model, pro- 
viding a distinction among the disk and the surrounding 
gas contributions. The corresponding geometrical factors 
are assumed t o be given by the isotropic interpretation of 
Netzer ( 1990? ) and the projection of a rotating disk, con- 
fined in a smaller region with respect to the other compo- 
nent. The inclination of the disk is estimated by compari- 
son of the BF symmetric distortion with that predicted by 
the reference model, as shown in Fig. 6, while Rblr is in- 
ferred by measurements of the optical continuum intensity, 
according to the empirical relationships fou nd in some re- 
cent reverberation rnappi ng campaigns (Kaspi et al. I . l200d 
200,4 iBentz et al. I ■ l20n6^ . 



A particularly important problem, involving the deter- 
mination of AGN physical properties from emission lines, 
resides in the line profile asymmetries. Several factors, 
such as partial obscuration, geometrical structure, or large 
scale non-virialized motions can produce asymmetric line 
profiles. Moreover, relativistic effects within the gravita- 
tional field of the SMBH give raise to asymmetries, espe- 



Figure 8: Degree of linear polarization at the radio frequency of 
1.4 GHz as a function of the inferred BLR inclination. The cross in 
the upper left corner gives the median uncertainty estimate. Polar- 
ization data are from the NVSS catalogue. 



dally in the high velocity wings of the profile, which are in- 
cluded in the calculations of the model by Ch en fc Halpern 
(198i). [n order to assess how much the asymmetric com- 
ponent affects our estimates of the velocity field, we intro- 
duced an asymmetry parameter: 



= /l3i/3(HWHMH/3), 



(16) 



expressing the relative contribution of the asymmetric com- 
ponent, with respect to the Gaussian component, in the 
profile of H/3 at its half-maximum level. As we show in 
Fig. 7, the asymmetric component gives a contribution 
to the FWHM which rarely exceeds the 10% level. The 
most extreme cases, where the asymmetric component be- 
comes larger than 20%, occur only in the range of very 
broad line emitting sources. Although this property is 
not particularly well represented in our sample, it echoes 
the observation of larger asymmetries in objects where 
FWHMh/3 > 4000 km s~^, which is arn ong the features 
identified by ISulentic et al. ] (l2000l . l200d) in their distinc- 
tion between Population A and B sources. Since there is 
evidence for different structural configurations in objects 
with broader and narrower emission lines, at least in the 
case of radio-loud sources, objects with the largest asym- 
metries are more problematic in their comparison with the 
reference model. 

5. 2. Implications 

Most of the results achieved in this work depend crit- 
ically on the choice of our reference model, which leads 
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Figure 9: The range of estimated geometrical factors for sources of 
different FWHMh/j. Open circles represent objects with symmetric 
profiles, while filled squares and triangles are for sources with posi- 
tive and negative asymmetries, respectively. The distinction among 
sources with different degree of asymmetry is required because ob- 
jects with large asymmetries are harder to compare with the theo- 
retical line profiles, predicted by the reference model. 



US to conclude that the BLR has a flattened component, 
that is commonly seen at i < 20° . In the case of radio-loud 
sources, nearly face-on disk structures are likely to produce 
a radio jet oriented along our line of sight and the resulting 
signal should be considerably polarized. Indeed, some ob- 
jects of our sample have been detected in the NRAO VLA 
Sky Survey (NVSS), which provides measurements of the 
radio flux and polariz ation at the frequency of 1.4 GHz 
( Condon et al. I . [l998f )F1 We compare the degree of linear 
polarization with our inclination estimates in Fig. 8. We 



find that a significant degree of linear polarization is de- 
tected in many objects and it appears to be an averagely 
decreasing function of i. 

A comparison of our mass determinations with the 
isotropic assumption allows us to study the properties of 
the geometrical factor within our sample. The situation 
depicted in Fig. 9 clearly indicates that significant effects, 
up to a factor ^ 40, should be expected and that they 
are more commonly observed in the range of sources with 
FWHMh/5 < 3000 - 4000km s-^ We find that the av- 
erage value of the geometrical factor for black hole mass 
determinations based on FWIIMh/3 is / = 6.51, not ver y 



different from the result achieved bv lOnken et al. I (|2004l ). 
who gave a geometrical factor / = 5.5 ± 1.9 for black 



^Polarization data are available at 
|http://www.cv.nr ao.edu/nvss/NVSSlist. shtmi] 



the web site 



Figure 10: SMBH accretion rates vs. FWHMh/3, using the same 
symbology as in Fig. 7. The cross in the upper right corner of the di- 
agram is the median uncertainty of measurements. We do not detect 
either very high accretion rates or systematic trends associated to the 
line profile widths. Our estimates do not appear to be dramatically 
affected by the presence of asymmetric profile components. 



hole mass estimates based on the emission line dispersions. 
We should note, ho wever, that a previous investigation by 
Collin et al. 1 (|2006^ applied composite kinematical models 
to the BLR and compared several estimates of the central 
engine mass, suggesting a distinction among the geomet- 
rical factors needed to correct the velocity fields based on 
the emission line dispersion and FWHM. 

Applying our modified scheme to the calculation of the 
central engine's dynamical properties, we find that a cor- 
rection for BLR inclination increases our estimate of the 
SMBH mass. A direct consequence of this effect is that 
we do not detect dramatically high accretion rates, with 
respect to the corresponding Eddington limits. In particu- 
lar, as it is illustrated in Fig. 10, although we do observe a 
trend where objects with narrower emission lines have typ- 
ically higher accretion rates, this property does not appear 
to hold in the form of a real anti-correlation, as it could 
look like, when the BLR geometrical properties were not 
properly considered. Comparing our estimates of accre- 
tion rate with the temperature parameters A, inferred by 
means of the BP analysis and plotted in Fig. 11, we find 
that the occurrence of high accretion rates is generally as- 
sociated to lower values of A (i. e. higher plasma temper- 
atures), but we detect no particular relationship with the 
effects that undermine the reliability of the BP method, 
suggesting that more physical processes should be involved 
in this problem. 
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Figure 11: The temperature parameter A compared with the es- 
timated accretion rate, in objects with different BP classification 
(symbols are the same as in Fig. 2). Lower values of A are observed 
with increasing Eddington ratios. 
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